Pattern Selection in a Phase Field Model for Directional Solidification 
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A symmetric phase field model is used to study wavelength selection in two dimensions. We study 
' the problem in a finite system using a two-pronged approach. First we construct an action and, 

minimizing this, we obtain the most probable configuration of the system, which we identify with the 
selected stationary state. The minimization is constrained by the stationary solutions of stochastic 
evolution equations and is done numerically. Secondly, additional support for this selected state is 
obtained from straightforward simulations of the dynamics from a variety of initial states. 
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The problem of pattern selection is to predict which pattern, out of a number of possible stationary patterns, will 
be selected under given experimental conditions. This problem has been of great interest for many years in many 
fields such as biology, chemistry, engineering and physics 0] where a large variety of reproducible patterns can be 
formed under various external conditions which drive the system away from thermal and mechanical equilibrium. 
Important examples of pattern selection are directional solidification and eutectic growth where the interface between 
the ordered and disordered phases can take a periodic cellular pattern. From the experimental point of view, the 
interface seems to have a well defined periodicity in both directed one phase solidification 0, HI m an d m directed 
eutectic solidification 0,0]. On the theoretical side, there is much disagreement with some researchers claiming that 
a unique wavelength is selected 0, H, IS E3- EH while others say that the wavelength of the final pattern is accidental 
fl2l flij . A lar ge s cale simulation on a noisy Swift Hohenberg equation 01 which models fluid convection near onset 
has been done |15| to study the selection of periodic patterns. For this system, a potential Tsh exists and minimizing 
this gives the unique selected stationary pattern for any initial state, exactly as for a system approaching equilibrium. 
' In this paper, we use a numerical and a theoretical approach to argue that in the presence of noise, there is a 
wavelength selection for the cellular pattern formation in the directional solidification problem, see Fig. 1. More 
precisely, we use a phase field model with additive stochastic noise developed for growth processes such as single 
phase directional solidification driven out of equilibrium by an external moving temperature gradient (lil UtI Il8| . 
This causes invasion of a disordered by an ordered phase with non potential dynamics. Since we are unable to solve 
the problem analytically, we use numerical computations which limit us to rather small system sizes. However, we 
believe that this is sufficient to demonstrate our conjecture of a unique selected state. In addition, we construct an 
action and, minimizing this, we obtain the most probable configuration of the system, which we identify with the 
selected stationary state. The minimization is constrained by the stationary solutions of stochastic evolution equations 
which are calculated numerically. 



PHASE-FIELD MODEL 



To carry out our simulations, we use a phase-field model for directional solidification based on two continuum fields 
to describe the phases of the system. These fields vary slowly in the bulk region and rapidly near the solid-liquid 
interface. The equilibrium system is described by a free energy functional T(c, ip) |l(| where 
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FIG. 1: A schematic representation of the cellular pattern in directional solidification. AT = (T — T m )/T m where T m is the 
melting temperature of the pure solid. 



Here, ip is a nonconserved order parameter field which has one value in the solid and another in the liquid and c 
is a dimensionless diffusion field proportional to the impurity concentration. AT = (T — T m )/T m where T m is the 
melting temperature of the pure solid when c = 0. Aip is the discontinuity in ip between solid and liquid, and K c 
are phenomcnological constants accounting for the free energy cost of spatial variations in ip and c. The Kobayashi 
function, h(ip) 19], is introduced for analytic and computational convenience to simplify the problem, despite the 
apparent additional complexity |20l El|. It has essentially the same effect as ip itself as h(ip) is chosen to have the 
properties (i) h(±l) = ±1 and (ii) dh/Sip — when ip = ±1. This choice ensures that the equilibrium values of ip 
are ip eq = +1 for AT > and ip eq = — 1 for AT < independently of the values of AT and c so that Aip = 2. The 
interface between the ordered AT < and the disordered AT > phases is the region over which the fields ip and c 
change from one bulk value to the other and has a width £ = 0(1), we define its position (x, z) by ip(x, z ) = 0- 

To impose the motion of the interface, we take AT(r, t) = AT(z — vt) where v = vi is the externally imposed 
pulling velocity. Stationary states can exist only in the frame moving with the interface when AT(r, t) — AT(z). 
In this frame, the simplest possible dynamics consistent with the macroscopic conservation law for c are Langevin 
equations 22] 
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where and r c are the mobilities of the fields ip and c respectively. Complications such as fluid flow are not 
considered. Note that the non dissipative terms vdipjdz and vdc/dz cannot be written as derivatives of a potential 
so that the dynamics of Eq. (|5J is non potential. This makes the construction of a differentiable Lyapunov functional 
for such situations difficult but such a procedure has been implemented successfully in a similar situation |23l l24l] . 
The stochastic noises r]^(r,t) and r] c (r,t) have probability distributions 
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When v = in Eq. the system approaches its stationary equilibrium state given by the minimum of T{c, ip) 
[i3,|2^]. G?(r, r') obeys V 2 G(r, r') = — S(r — r') with appropriate boundary conditions and G _1 (r, r') = — V 2 (5(r — r') 
so 

(^(r,t)77^(r',0> = 2T^5(r - r')6(t - t') 

( Vc (r,t) Vc (r',t')) = 2T' c G- 1 (ry)S(t-t'). (4) 
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The form of Eq. Q is required by the dynamics of the nonconserved field ip and the conserved field c [22]. If 
c = r^, jC , then riTp and r\ c obey a fluctuation dissipation theorem but this is not necessary and makes no quali- 
tative difference. In our numerical simulations, the noise strengths V, c are determined by the widths of the noise 
distributions. The exact form of the Green's function G(r,r') in Eq. ||3j) is needed for numerical computations. 



PATH INTEGRAL AND STATE SELECTION 



The model defined by Eqs. @ and © describes the equilibrium of the two phase system and the invasion of one 
phase by another in as simple a way as possible. Consequently, we take the noise strengths and mobilities T c ^ to be 
field independent constants. This simplified model does describe equilibrium and near equilibrium dynamics correctly 
and is a good phenomenological model far from equilibrium. In our model, the temperature gradient AT(r, t) is taken 
to not affect the noise strengths. This may not be entirely realistic but is sufficiently simple to permit analysis of the 
model and does describe the essential features of directional solidification. The great advantage of our model is that 
complications such as Landauer's "blowtorch theorem" j2f| are not applicable despite the presence of a "temperature" 
gradient AT as the noise strength is taken to be independent of T. 

We formulate the problem in terms of an action S(tp, c) whose minimum gives the most probable or selected state. 
The joint probability distribution V(4>, c) is 



P(ip, c) — J Vr) 1 pVi] c 'P{r)ip)'P{r) c )S(ri^ - tJj + Q^,)8{i] c - c + Q c 
= exp —S(ip, c) 

using Eqs. (|2I3|> . The configuration (-0, c) minimizing S(ip, c) is the most probable state at least for weak noise 1261127 
The explicit expression for the action is 
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where, Gij> and Q c are defined in Eq. J2J) and G(r, r')is the Green's function. Note that this quadratic form for S(ip, c) 
is very similar to the actions proposed previously [ill I23I Eil |2t| . The action of Eq. © can be obtained from the 
Martin-Siggia-Rose (MSR) formalism |3(J, |3l| by integrating out the auxiliary fields 'ip and c. Since the MSR action 
contains all physical response and correlation functions, including those at equilibrium. The resulting action of Eq. 
© is also a complete description of the dynamics and of an equilibrium system when = —8T /Sip and similarly 
for Q c which ensure that the dynamics of Eq. J5| is satisfied. We are interested in stationary patterns at very large 
times t. A reasonable assumption is to argue that the pattern becomes very close to the selected stationary pattern 
at some large but finite to and remains in this stationary configuration ^(r); c(r) for t > to when tp — = c. Then, 
the action is dominated by the infinitely long time interval T = t — to — > 00 and becomes 
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the above Eq. is our essential result which should be regarded as a conjecture as a real proof eludes us. 

As discussed below Eq. JSJ, the configuration ip(r), c(r) which minimizes the action S(ip,c) in Eq. J7J maximizes 
the probability V(ip, c )- This is the selected stationary state if dynamical paths exist between any intermediate state 
and the selected state in analogy with the ergodic theorem. This is almost certainly true but is very hard to verify by 
simulations as metastable states in the Eckhaus band [32] ] have very long lifetimes with weak noise even for our very 
small systems. This is similar to metastable states separated from the equilibrium state by large free energy barriers. 
Such effects lead to apparent history dependence of the final state in experiments |3^| and in theoretical investigations 
without noise [34| ■ The essential difference in the case studied here is that the dynamics of Eq. (J2J) is non potential 
while the approach to equilibrium is governed by potential dynamics A naive minimization of the action S(ip, c) 
of Eq. JHJ yields the deterministic evolution equations which are incorrect with noise. The minimization must be done 
subject to the constraints that ^( r ;^) an d c(r, t) are solutions of Eq. © which we perform numerically. In principle, 
this can be done analytically in the MSR formalism [3lJ, as used in studies of the KPZ system [35| . 
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FIG. 2: Evolution of cellular pattern from an initial random interface to final state qf = 0.02 with velocity v = 0.20 in absence 
of noise, and L x — 500. 

NUMERICAL PROCEDURE 

To test our theory numerically, we generate from simulations of Eq. J2J all quasi-stationary patterns ^>(r), c(r) with 
noise and finite v and compute 4T'S(ip,c)/T from Eq. JJJ. In the simulations, parameters similar to those describing 
experiments on the liquid crystal system 4-n-octylcyanobiphcnyl (8CB) [4],|36J are used: r c =2r^, = 1, 7 = 0.73, 
K 4 , = 1.5, K c = and the pullingspeed v = 0.20. We take AT{z) = -T for z < -W, +T„ for z > W and Gz 
for |*| < VF with G = T /W [lfil Til Haj| . The temperature T = 0.38 so that AT(z) is greater than the liquidus 
temperature for z > W and below the solidus for z < —W. The simulations are carried out with the width of the 
temperature variation W = 50, the size of the simulation box L z = 150, lattice spacing Ax = 1.0 and timestep 
At = 0.10. The system widths are 250 < L x < 1000 with periodic boundary conditions %jj(x + L x ,z,t) = ip(x,z,t) 
and c(x + L x , z, t) = c(x, z, t) which restricts the allowed values of the q vector of a periodic variation in the interface 
between the ordered and disordered phases to be q x = n/L x with n = 1,2, ■ ■ ■ L x . The Green's function G(r,r') 
satisfies G(x, 2; x' , z') = G(x + L Xl z; x', z') and G(x, 0; x' , z') = = G(x, L z + 1; x', z') and is obtained by numerically 
inverting G(p), the discrete Fourier transform of G(r, r'). This is used to calculate AS = 4T'(S n — S m )/T where the 
integers n,m — q x L x label the patterns. 

RESULTS AND DISCUSSION 

We perform many simulations of Eq. (J2J) from various initial conditions in order to verify if a single pattern, out of 
a number of possible stationary patterns, is indeed selected in the presence of weak noise. 

First, simulations are done in the absence of noise 77^, = — r\ c but with initial values of ip = e in such a way to form 
an initial random interface from a uniformly distributed random variable and c — for convenience. For system widths 
L x = 250, 500 and 1000 the system always eventually reaches a periodic stationary state with qf = n/L x = 0.020, 
where n is the number of cells. For L x = 420, the final periodicity is qf = 8/420 « 0.01905. These results imply that 
the selected wavevector is the nearest allowed qf to the preferred value for L x — > 00. In Fig. we show the time 
evolution of the random interface for L x — 500, after t = 10 5 time steps the system reaches a pattern with 10 cells, 
i.e, q f = 0.02. 

Second, we run simulations with an initial periodic interface of the form £(x) = (ocos(27r<7o20 with various wavevec- 
tors qo = no/L x and amplitude Co similar to that of the final state of the previous simulation, again in the absence 
of noise. The initial interface is located at the center of the simulation box of size L x x L z . Simulations with various 
allowed values of qo and L x = 500 yield final stationary periodic states with qf = 0.020, 0.022 and 0.024 which 
determine the Eckhaus stable band of wavevectors for our parameter values. In an infinite system, the Eckhaus band 
has a continuum of q values which is reduced to a small finite number by the finite width L x — 500. When the size is 
doubled to L x = 1000, we find 7 wavevectors in the Eckhaus stable band from q x — 0.019 to 0.025 as expected. 

The next simulations study the effect of noise on the evolution of an initial periodic state with wavevector qg = 0.016 
outside the Eckhaus stable band. Without noise, the evolution is by tip splitting of alternate cells to a final wavevector 
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FIG. 3: Evolution of cellular pattern of initial periodicity qo — 0.016 to final state qj = 0.024 with velocity v = 0.20 in absence 
of noise, and L x — 500. 
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FIG. 4: Time evolution of same initial pattern as Fig. Q with stochastic noise of strength rj = 0.10 to final state qf = 0.020. 

of qf — 0.024, as shown in Fig. (J3J. This configuration of 12 cells is different from the expected wavelength 
of qf — 0.020 as obtained from random initial conditions. With applied noise rj^(r,t) independently uniformly 
distributed with —77 < tj^(t, t) < +rj with maximum noise amplitude 77 = 0.10, the evolution is by a totally different 
route to the expected state with q$ = 0.020 as shown in Fig. (@J. This indicates that noise is essential in the selection 
of the final state and that only one wavenumber in the Eckhaus band is truly stable in the presence of stochastic noise. 
We have also performed many simulations, not shown here because of space limitations, with several different initial 
configurations all of which eventually go to the same state with qf — 0.020 with only a single exception discussed 
below. 

Now that we have very strong evidence for selection of a unique wavelength from almost all initial patterns except 
for some qo in the Eckhaus band, we perform a simulation with an initial pattern = ^o co ^(^Qo x ) with qo = 0.024 
with noise to study the stability of modes in the Eckhaus stable band [l|,|32j. In the absence of noise, all modes with 
qo in this band are stable. We impose an interface £(x), evolve it in absence of noise for 10 4 time steps to allow higher 
harmonics to develop and then evolve the system with noise of strength rj = 0.20 for 2 x 10 5 time steps with results 
shown in Fig. ©. The interface evolves by the gradual loss of one cell to a final configuration of 11 cells or qf = 0.022 
which remains stable for at least 5 x 10 7 steps. When a brief burst of very strong noise with 77 > 0.25 is applied and 
the system then allowed to evolve under weak noise, the expected qs — 0.020 state results. This is similar to the 
history dependent effects observed |3^| and discussed theoretically in the absence of noise [34| • 

In order to verify our conjecture, we calculate the action from Eq. Q of the three possible states in the Eckhaus band. 
As shown in Fig.©, the action difference AS(ip,c) oc rj 2 so that, for the deterministic case, AS = 0. A finite AS is 
due to finite noise implying that the selection of a particular pattern is due to noise driven fluctuations. Our L x = 500 
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FIG. 5: Evolution of an interface with go = 0.024 with noise of strength r\ — 0.10 to a final state qj — 0.022. 



system has only three quasi-stationary periodic states in the Eckhaus stable band [32j with q x = 0.020, 0.022, 0.024 
corresponding to n = 10, 11, 12 cells. The qs — 0.020 minimum action state is conjectured to be the unique selected 
state. An outstanding question is whether pattern selection also holds in the thermodynamic limit L x — oo when the 
set of stationary states in the Eckhaus stable band is continuous. Comparing with experiment requires extending the 
calculation to 3D to study the selected stationary structures of the 2D interface |l6L l37| . 
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FIG. 6: Action differences AS = 4r'(<Sn - S w )/T and 4r'(5i 2 - S w )/T from Eq. in a system width L x = 500 where the 
integers n = q x L x in <S n . This shows that the pattern with wavelength qs = 0.020 has minimum action. The noise is uniformly 
distributed with —rj < rj(r, t) < +rj. 



These numerical simulations are consistent with the hypothesis of a unique selected stationary state with a definite 
periodicity q$ = 0.020. This is the noise induced stationary state which is reached from all initial states except from 
the periodic state with q = 0.022 in the Eckhaus stable band. This single state is stable with weak noise up to 
our longest simulation time. If a unique selected state exists, then only one of these should be stable to weak noise 
but some are stable within our computation time and it is impractical to perform sufficiently long runs to detect an 
instability toward the hypothetical minimum action selected state with qs = 0.020. This study finds that all states 
with this single exception do evolve to the minimum action state and, in addition, provides some insight into the 
relative stability of states in the Eckhaus band. On balance, there is overwhelming evidence for our hypothesis of a 
unique selected state (tp, c) determined by the minimum of the action iS(^, c) of Eq. lf?ll. 

One important caveat must be made. In analogy with a finite equilibrium system |38(, we expect all states in our 
finite system to have finite lifetimes as the stochastic noise causes transitions between them. We expect that true 
stationary states exist only in the thermodynamic limit of L — > oo when the lifetime of the selected state should be 
infinite but this limit is extremely difficult to approach numerically for well known and obvious reasons. If a finite 
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size scaling theory existed for these dynamical processes, the L — > oo limit could be studied by the methods of this 
letter. However, to the best of our knowledge, such a scaling theory does not exist. 
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